Coherent adiabatic theory of two-electron quantum dot molecules in external spin baths 
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We derive an accurate molecular orbital based expression for the coherent time evolution of a two-electron 
wave function in a quantum dot molecule where the electrons interact with each other, with external time depen- 
dent electromagnetic fields and with a surrounding nuclear spin reservoir. The theory allows for direct numerical 
modeling of the decoherence in quantum dots due to hyperfine interactions. Calculations result in good agree- 
ment with recent singlet-triplet dephasing experiments by Laird et. a/.[Phys. Rev. Lett. 97, 056801 (2006)], 
as well as analytical model calculations. Furthermore, it is shown that using a much faster electric switch than 
applied in these experiments will transfer the initial state to excited states where the hyperfine singlet-triplet 
mixing is negligible. 

PACS numbers: 73.21.La, 78.67.-n, 85.35.Be, 78.20.Bh 



It is now well recognized that the hyperfine interaction is 
one of the main sources of decoherence in few-electron quan- 
tum dots. This interaction, originall^considered in metals by 
Overhauser more than 50 years ago couples the electronic 
spin states through weak nuclear spin interactions with an or- 
der of ~ 10® surrounding nuclei yj]. Recently, this coupling 
has received considerable interest through the demonstration 
of controlled single electron manipulation 0] which opens 
for real quantum information processing based on electronic 
spin states in quantum dots jsf]. Procedures to minimize 
or control the hyperfine interaction is therefore vital for the 
functioning of any quantum dot based information processing 
technology. Such mechanisms also contain novel aspects of 
spin de- and re-phasing of quantum systems interacting with 
a large spin bath. 

Recently it was demonstrated in experiments |@, 0] with 
two-electron quantum dot molecules that the magnitude of 
the hyperfine interaction is consistent with a random magnetic 
field strength of a few mT. The nuclear field-induced singlet- 
triplet coupling leads to a spin dephasing of an initially pre- 
pared singlet state within 1 — 10 ns. These experiments utilize 
fast adiabatic electric switching techniques which transform 
the ground state from a two-electron ionic state in one dot to 
a covalent state with one electron in each dot. 

The experiments have been analyzed in detail theoretically 
based on various model Hamiltonians For small tunnel- 

ing coupling the effective two-electron Hilbert space amounts 
to the four possible covalent spin states (a singlet and three 
triplet states). This may be further reduced to two states by 
exposing the molecule to a magnetic field which decouple the 
two states with nonzero magnetic quantum numbers. Within 
this approximation it was shown that the hyperfine interaction 
induce a spin saturation which saturates sensitively as a func- 
tion of the exchange coupling and the hyperfine coupling . 
Th ese p redictions were confirmed experimentally by Laird et. 
al. IM- 

In the present Letter we develop a full coherent model of 
the two-electron spin dynamics which includes the hyperfine 
interaction on equal footing with the time dependent external 
fields. The results may be directly compared with the experi- 
mental results. The theory not only validates the effective two- 



level models in the presence of an external magnetic field, but 
also demonstrates predictive saturation values in the absence 
of an external magnetic field. Furthermore, it will be shown 
that decreasing the switching time an order of magnitude may 
lead to controlled diabatic transfer ifTHl to excited states where 
the singlet-triplet mixing can be neglected. 

Our starting point is the Hamiltonian of two interacting 
electrons in a double quantum dot with dot separation d as 
recently applied in studies of electron structure as well as in 
studies of photon induced controlled transport. lfl2l[l3i[l4ll . 



H = h{ri) + h{Y2) + 
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Here ri 2 are single-particle coordinates, ^ is an electric time 
dependent field applied along the inter-dot axis and B^xt is an 
external magnetic field perpendicular to the dot. The material 
parameters are taken as those of GaAs, with effective mass 
m* — 0.067me and relative permittivity — 12.4. The gy- 
romagnetic ratio for GaAs is 7e = <7* with g* ~ —0.44. 
The confinement strength is set to = 1 meV and the inter- 
dot separation to d ~ 130 nm, which are realistic experimen- 
tal values Rfisll. 

Fig. [U shows the energy spectrum obtained from diagonal- 
ization with Bpxt = and Bpxt = 200 mT (inset). The 
spectrum of the system has been explained in detail else- 
where ifisll . For Bext = symmetry about the y-axis is con- 
served and we show only the singlet and triplet (iris = 0) 
eigenstates corresponding to the three lowest energy bands 
with even reflection symmetry. This symmetry is broken when 
the external magnetic field is applied and all states are shown. 
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When Bext 7^ the magnetic sub-levels split (not shown). In 
addition to the splitting of the spin states (anomalous Zeeman 
effect) the spectrum is also changed by the Zeeman (L^) and 
diamagnetic (-Bg^-j) terms, as seen clearly in the inset. The 
most important effects of these terms are the modified singlet- 
triplet splitting J, the modified anti-crossing energy difference 
and the splitting of the second band according to the sign of 
the angular momentum expectation value, {L^)- 

It is particularly worth noting that the energy spectrum ex- 
hibit several near degenerate anti-crossing regions where the 
coupling strength can be tuned in experiments through ad- 
justable gate voltages and switching times |10!]. Restricted 
by conservation of total spin, the states can couple dynami- 
cally as the electric field varies. The relative coupling strength 
from our model is shown in the lower panel of Fig. [T] The 
strongest coupling strength is seen between the singlet ground 
state and first excited states at —0.013 mV/nm. This puts a 
limit on the switching time through the region enclosed by a 
green circle for adiabatic time development along the initial 
singlet ground state, marked as T' in Fig. [T] On the other 
hand, a very rapid transfer can lead to diabatic development, 
to be discussed later. 

The few and well defined states resulting from the 
present diagonalization suggest that the time evolution is 
most precisely described in an adiabatic basis of instant 
eigenstates of the time-dependent Hamiltonian of Eq. ([l]), 
H{t)x{r i,r2;S.) = e(0x(ri,r2;C), where the energies e(f) 
depend parametrically on the time-dependent electric field. 
We expand the wavefunction in these basis states. 
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where k runs over all basis states. Projecting onto each basis 
state, the following expression for the time evolution of the 
amplitudes is obtained. 



{Xk\X\xj) 



Cjit)+iek{Ockit), (4) 



with X = xi + X2- 

Additional terms can readily be included as extra matrix 
elements in the expression above. To study hyperfine inter- 
actions, spin couplings for ~ 10^ nuclear spins surrounding 
the electrons must be included. These evolve in time, but on a 
much longer time scale than we will consider here. We there- 
fore use the quasistatic approximation (Hld^, neglecting their 
time dependence. In addition, the large number of spins justi- 
fies a semiclassical approximation [17], where all the nuclear 
spins are described by a single classical magnetic field. The 
hyperfine interaction is then given by 
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where 5, is the spin operator of electron i. Gener- 
ally the direction of nuclear magnetic field Bn is ran- 
dom (no polarization) and the magnitude varies accord- 
ing to a normal distribution about zero, P{Bm) — 
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FIG. 1: (color online) Upper panel shows a few of the lowest en- 
ergy levels as a function of electric field in the a;-direction, corre- 
sponding to y-symmetric states. Inset shows the effect of an exter- 
nal magnetic field (200 mX) on the spectrum. Lower panel shows 
the relative coupling strengths of selected states; Full blue curve is 
(1, S|X|0, S) between the two lowest singlet states (full blue curve), 
full red curve is (1, To |X|0, To) between the two lowest triplet states 
and (3, 5*1X11, S) dashed blue curve is between the second lowest 
and third lowest singlet. Note that the rris = ±1 triplet states are not 
shown. 



1/(2^B2^,J2 exp {~Bn ■ Bn/2BI^^) M- 5„uc can be de- 
termined by experiments and is of the order of 1 mT |3l . The 
precise spatial variation of the nuclear magnetic field Bn is in 
general unknown, but also of less importance. The essential 
feature in the spin dephasing mechanism is the difference in 
effective magnetic fields between the two dots. The simplest 
way to represent this is by a step function. 
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N = 




, for a; > 
, otherwise 



(6) 



This induces coupling between the singlet and triplet states 
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FIG. 2: (color online) Upper panel: Averaged time evolution of sin- 
glet component with Bext ~ 200 mT. Bn/J ~ 0.4 (full blue line) 
and Bn/J = 0.91 (full red line, downshifted 5%). Squares and 
circles are experimental data taken from fioll . Prediction from two- 
level model shown as dashed lines. Semitransparent dots are actual 
numerical data, full lines are obtained from a smoothed spline inter- 
polation. Lower panel: Bm/J = 0.91, but without external mag- 
netic field. Prediction from four-level model (dashed line), scaled to 
match at t = 0. 



and between the different triplet states. 

From Eqs. ( I1I3I6I I the time evolution of the wavefunction at 
zero electric field, restricted to the four lowest energy states 

{|5),|To),|T_),|r+)}, becomes. 
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This expression is identical to four-level models previously 
considered |@]. We have here, however, excluded the inter- 
triplet couplings as this will allow us to obtain an analyti- 
cal solution. When the external magnetic field is sufficiently 
strong, the rus = ±1 triplet components effectively decouple, 
and we are left with a two-level system defined by the upper 
left part of the four-level matrix. Furthermore, since the triplet 
sub-levels are degenerate, the four-level matrix for B^^t = 
may be represented by an effective "radial" two-level model. 
In both cases we obtain the time evolution of the singlet coef- 
ficient 



|c5(t)P = l- 



4^2 



4B2 + J2 
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J2 



(8) 



with B = B, for B^xt > J and S = ^ Bl + B^ + B^ for 

Bext — 0. 

In Fig. 12] we show our results together with experimen- 
tal results of Laird et. al ifioll using Bext — 200 mT (up- 



per panel). To calculate the singlet-triplet dephasing we 
start out in the singlet ground state at large electric field, 
—0.024 mV/nm. The field is then switched adiabatically to 
zero within 1 ns and kept at zero a variable period of time, 
ts up to 16 J/h ns, before being switched back to its original 
value. The procedure is repeated a number of times, with a 
random nuclear magnetic field drawn from a Gaussian distri- 
bution. Finally, the average of the singlet correlator is com- 
puted, Ps = ^/NJ2i |cs(is)P- The full fines in the upper 
{Bext = 200 mT) and lower panel (B^xt = 0) of Fig. |2] 
are produced from a sample size of = 1000 different nu- 
clear fields. For clarity, we have interpolated the numerical 
data using splines with a smoothing requirement. The actual 
numerical points are shown as semi-transparent dots. Also 
shown, as dashed lines, are predictions from the two-level the- 
ory. An excellent agreement between the theoretical results 
are noted which indicate that effects of the electrical switch, 
excited states and geometry of the potential are of less im- 
portance in this case. The present results are also compared 
with experimental data, shown as dots and circles in the up- 
per panel. We also observe a very good agreement with the 
experimental results. Here we have varied i?nuc as opposed 
to J which may be varied in experiments by tuning the gate 
voltages. As verified by the two-level models the dephas- 
ing process mainly relies on their ratio. It should be noted 
that results in the upper panel have been scaled according to 
Ps{t) = 1 - V"(l - P^(t)), where V = 0.40 is a visibil- 
ity parameter determined in the experiments, Ps {t) is the ex- 
perimental averaged correlator while Pg (t) is the theoretical 
averaged correlator ifioll . 

Setting the external magnetic field to zero leads to increased 
dephasing, since the = ±1 states are now coupled to the 
initial singlet state. This is indeed what we observe in our 
simulations, as shown in the lower panel of Fig.|2l The four- 
level (effective two-level) model yields the result shown as a 
dashed line, scaled to match the numerical data at ts = 0. At 
large evolve times, these are in good agreement. Numerical 
solution of the full four-level matrix suggests that the slight 
discrepancy around tg = 2 is due to the neglected inter-triplet 
couplings in Eq. (|7]i. In contrast to the two-level case, we 
observe in the four level case a ~10% dephasing occurring 
during the switch or more precisely between the two passages 
of the avoided crossing circled out in Fig.[T] This is a result of 
a much more involved dynamical interplay between the states 
and result in the present theory to a reduction of Pg at t,, = 0. 

A completely new feature can be studied by introducing an 
ultrafast switching function which leads to diabatic passage 
through the anti-crossing, highlighted by the green circle in 
Fig.[T| This will transfer the system to the first excited singlet 
state where at zero electric field the singlet-triplet splitting 
is approximately 100 times greater compared to the ground 
state. Important to note is that states in the second energy 
band are like the ground state covalent states lfl3ll with one 
electron in each dot. This makes the states well suited for 
single-electron gate operations unlike the states in the third 
energy band which resemble ionic states. In Fig.[3]the evolu- 
tion of the three lowest energy singlet states are displayed as 
the system is rapidly switched from positions marked I to II 
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FIG. 3: (color online) The time evolution of the three most promi- 
nent states during the fast switch proceedure. About 95 % of the total 
norm is represented by these states, which are 15", 0) (blue dashed), 
15*, 1) (red full) and |S,2) (black dash-dotted). B^uc = 1 mT 
([1,1,1]). Bext ~ 0. Left panel shows evolution during rapid switch 
(1 ps), center panel shows adiabatic passage to first excited state and 
evolution during dephasing period. Right panel shows rapid switch- 
ing back to the one-well configuration (1 ps). Singlet-triplet coupling 
is weak for the excited singlet state, and thus only ~ 5 % is lost 
(mainly to other singlet states during switching). 



in Fig.[Tl using a nuclear magnetic field of 1 mT, and zero ex- 
ternal magnetic field. From there it is transferred adiabatically 
from II to lib. The leftmost panel shows the transfer of popula- 



tion from the ground to first excited state during a 1 ps switch, 
with around 10% of the population being further transported 
to the second excited state. In the middle panel, the system is 
switched adiabatically to zero electric field during 2 ns, and 
left to evolve here for 50 ns. The singlet-triplet coupling is 
completely suppressed, and when the system is switched back 
to position I, 95% of the initial singlet population is regained, 
the rest having vanished mostly to higher excited states dur- 
ing passage through the anti-crossings. The simulation was 
repeated a number of times with increasing nuclear magnetic 
field strength up to 10 mT, attributing only negligible changes 
to the dynamics. We note that by applying optimal control 
schemes the transition to excited states n iay be achieved with 
near 100% transition probabiUty ifTH [ll [isll . 

In conclusion, we have applied a first principle molecular 
orbital framework to accurately calculate the time develop- 
ment of the electronic states of a two-electron quantum dot 
molecule. For the first time the hyperfine interaction has been 
taken into account on the same coherent level as the time vari- 
ation of the external electromagnetic fields. Calculations have 
displayed a very good agreement with experiments and a pre- 
viously developed two-level model for experiments performed 
in static external magnetic fields. We have also performed full 
numerical calculations in absence of external magnetic field 
and predicted the dephasing dynamics by an effective four 
level model. Finally, we have pointed towards a new solu- 
tion to the spin-dephasing problem by controlling and apply- 
ing transitions between excited states. 
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